home *** CD-ROM | disk | FTP | other *** search
/ Mac Magazin/MacEasy 2 / Mac Magazin and MacEasy Magazine CD - Issue 02.iso / Themen Mac Magazin / Multimedia / Grafik-&QT Tools / NIH Image 1.53 (non-fpu) / Macros / Input⁄Output Macros < prev    next >
Text File  |  1994-01-04  |  13KB  |  524 lines

  1. var
  2.   CharLoc:integer; {Global variable, initially zero}
  3.  
  4. macro 'Save using Time as Name…';
  5. {Note: Colons are not allowed in file names.}
  6. var
  7.   year,month,day,hour,minute,second,DayOfWeek:integer;
  8. begin
  9.   GetTime(year,month,day,hour,minute,second,DayOfWeek);
  10.   SaveAs(year-1900:2,'-',month:2,'-',day:2,
  11.          '/',hour:2,'-'minute:2,'-',second:2);
  12. end;
  13.  
  14.  
  15. macro 'Open with selection… [O]';
  16. begin
  17.   if nPics>0 then KillRoi; {Save Selection}
  18.   Open('');                {Prompt for file name}
  19.   RestoreROI;              {Transfer selection to new window}
  20. end;
  21.  
  22.  
  23. macro 'Save All…';
  24. {
  25. Saves all currently open images in a folder using '001', '002', etc.
  26. as the file names. The save file dialog box will be displayed once
  27. (and only once) so that you can specify the folder to save the files in.
  28. Leave the file name blank(e.g. SaveAs('')) to get a dialog box for each file.
  29. }
  30. var
  31.   n:integer;
  32. begin
  33.   RequiresVersion(1.45);
  34.   for n:=1 to nPics do begin
  35.     SelectPic(n);
  36.     SaveAs(n:3);
  37.     {Export(n:3);}
  38.   end;
  39. end;
  40.  
  41.  
  42. macro 'Import FITS…';
  43. {
  44. Imports 8 and 16-bit FITS images. Refer to "FITS:A Flexible Image
  45. Transport System", Astronomy and Astrophysics Supplement
  46. Series 44, 1981, 363-370.
  47. }
  48. var
  49.   width,height,offset,HdrPid:integer;
  50.   recsize,maxrecs,maxlines,loc,i,line:integer;
  51.   BitPix,str:string;
  52. begin
  53.   RequiresVersion(1.52);
  54.   recsize:=2880;
  55.   maxrecs:=10; {Must be increased to handle headers longer than 360 lines.}
  56.   maxlines:=(recsize*maxrecs)/80;
  57.   width:=recsize; 
  58.   height:=maxrecs;
  59.   offset:=0;
  60.   SetImport('8-bit'); 
  61.   SetCustom(width,height,offset);
  62.   Import(''); {Read in header as an image, prompting for the file name.}
  63.   HdrPid:=PidNumber;
  64.   BitPix:=Concat(chr(GetPixel(108,0)),chr(GetPixel(109,0)));
  65.   if (BitPix<>' 8') and (BitPix<>'16') then begin
  66.     PutMessage('This macro only reads 8 and 16-bit FITS files');
  67.     SelectPic(HdrPid);
  68.     Dispose;
  69.     exit;
  70.   end;
  71.   str:=concat(chr(GetPixel(266,0)), chr(GetPixel(267,0)),
  72.     chr(GetPixel(268,0)), chr(GetPixel(269,0)));
  73.   width:=StringToNum(str);
  74.   str:=concat(chr(GetPixel(346,0)), chr(GetPixel(347,0)),
  75.     chr(GetPixel(348,0)), chr(GetPixel(349,0)));
  76.   height:=StringToNum(str);
  77.   loc:=0;
  78.   line:=0;
  79.   repeat
  80.     str:='';
  81.     for i:=1 to 3 do begin
  82.        str:=concat(str,chr(GetPixel(loc mod recsize, loc div recsize)));
  83.        loc:=loc+1;
  84.     end;
  85.     loc:=loc+77; {Skip to end of line}
  86.     line:=line+1;
  87.   until (line>=maxlines) or (str='END');
  88.   SelectPic(HdrPid);
  89.   Dispose;
  90.   offset:=recsize+recsize*((line*80) div recsize);
  91.   if BitPix=' 8'
  92.     then SetImport('8-bit')
  93.     else SetImport('16-bit Signed; Calibrate; Autoscale');
  94.   SetCustom(width,height,offset);
  95.   Import('');  {No prompt this time; Import remembers the name.}
  96.   if BitPix=' 8' then Invert;
  97.   FlipVertical;
  98. end;
  99.  
  100. macro 'Display FITS Header…';
  101. var
  102.   i,j,loc:integer;
  103.   HdrPid,RecordSize,maxlines,line:integer;
  104.   MaxRecords:integer;
  105.   wname,str:string;
  106. begin
  107.   RequiresVersion(1.53);
  108.   RecordSize:=2880;
  109.   MaxRecords:=10;
  110.   maxlines:=(RecordSize*MaxRecords)/80;
  111.   SetImport('8-bit'); 
  112.   SetCustom(RecordSize,10,0);
  113.   Import(''); {Read in header as an image, prompting for the file name.}
  114.   HdrPid:=PidNumber;
  115.   SetFont('Monaco');
  116.   SetFontSize(9);
  117.   wname:=concat(WindowTitle,' Header');
  118.   NewTextWindow(wname,500,500);
  119.   loc:=0;
  120.   line:=1;
  121.   repeat
  122.     ChoosePic(HdrPid);
  123.     str:='';
  124.     for j:=1 to 80 do begin
  125.        str:=concat(str,chr(GetPixel(loc mod RecordSize, loc div RecordSize)));
  126.        loc:=loc+1;
  127.     end;
  128.     SelectWindow(wname);
  129.     writeln(str);
  130.     line:=line+1;
  131.   until (line>maxlines) or (pos('END',str)=1);
  132.   SelectPic(HdrPid);
  133.   Dispose;
  134.  end;
  135.  
  136.  
  137. macro 'Import Image TIFF File…';
  138. {
  139. As an example of how to import a foreign file format, this macro reads
  140. the TIFF files created by Image. The format of an Image TIFF file
  141. is described in Appendix E of the Image manual.
  142. }  
  143. var
  144.   width,height,offset:integer;
  145. begin
  146.   width:=768; 
  147.   height:=1;
  148.   offset:=0;
  149.   SetImport('8-bit'); 
  150.   SetCustom(width,height,offset);
  151.   Import(''); {Read in header as an image, prompting for the file name.}
  152.   if not ((GetPixel(0,0)=77) and (GetPixel(0,0)=77)) then begin  {'MM'}
  153.     PutMessage('This is not a TIFF file.');
  154.     SelectPic(nPics); Dispose;
  155.     exit;
  156.   end;
  157.   width := (GetPixel(30,0)*256) + GetPixel(31,0);
  158.   height := (GetPixel(42,0)*256) + GetPixel(43,0);
  159.   SelectPic(nPics);  {The ID of the last window opened is equal to nPics.}
  160.   Dispose;
  161.   offset:=768;
  162.   SetCustom(width,height,offset);
  163.   Import('');  {No prompt this time; Import remembers the name.}
  164. end;
  165.  
  166.  
  167. macro 'Import Multiple Images per File…';
  168. {
  169. Imports a series of 256x256 images contained in a single file, in this
  170. case an NIH Image stack with an arbitrary number of 256x256 slices.
  171. }
  172. var
  173.   offset,i,PicSize,HdrSize,width,height:integer;
  174. begin
  175.   HdrSize:= 768;
  176.   width:= 256;
  177.   height:=256;
  178.   PicSize:=width*height;
  179.   offset:=HdrSize;
  180.   SetImport('8-bit');
  181.   for I:=1 to 100 do begin  {Macro will terminate at eof}
  182.     SetCustom(width,height,offset);
  183.     Import('');
  184.     offset:=offset+PicSize;
  185.   end;
  186. end;
  187.  
  188.  
  189. macro 'Import PET…';
  190. var
  191.   offset,i,PicSize,HdrSize,width,height:integer;
  192. begin
  193.   HdrSize:= 0;
  194.   width:= 128;
  195.   height:=128;
  196.   PicSize:=width*height;
  197.   offset:=HdrSize;
  198.   SetImport('8-bit');
  199.   for I:=1 to 100 do begin  {Macro will terminate at eof}
  200.     SetCustom(width,height,offset);
  201.     Import('');
  202.     offset:=offset+PicSize;
  203.   end;
  204. end;
  205.  
  206.  
  207. macro 'Convert Files…';
  208. {
  209. Converts a set of raw data files(all in the same folder) with names
  210. in the form raw.001, raw.002, etc to TIFF or PICT.  As long as the
  211. converted files are saved in the same folder, you should
  212. only see two file dialog boxes(one for the first Import and one for
  213. the first SaveAs).
  214. }
  215. Var
  216.   i,nFiles:integer;
  217. begin
  218.   nFiles:=GetNumber('Number of files:',5);
  219.   for i:=1 to nFiles do begin
  220.     Import('raw.',i:3);
  221.     SetPicName('file',i:3);
  222.     SaveAs;
  223.     Dispose;
  224.   end;
  225. end;
  226.  
  227.  
  228. macro 'Import IPLab File';
  229. var
  230.    width,height,offset:integer;
  231. begin
  232.    width:=100; 
  233.    height:=1;
  234.    offset:=0;
  235.    SetImport('8-bit'); 
  236.    SetCustom(width,height,offset);
  237.    Import(''); {Read in header as an image, prompting for file name.}
  238.    width := (GetPixel(8,0)*256) + GetPixel(9,0);
  239.    height := (GetPixel(12,0)*256) + GetPixel(13,0);
  240.    Dispose;  
  241.    offset:=2120;  {The IPLab offset}
  242.    SetImport('16-bit Signed; Calibrate; Autoscale');
  243.    SetCustom(width,height,offset);
  244.    Import('');  {No prompt this time; Import remembers the name.}
  245. end;
  246.  
  247.  
  248. procedure ShowBioRadInfo;
  249. {Displays the contents of the 480(?) byte header at}
  250. {the end of Biorad MRC 600 Z Series files.}
  251. var
  252.   MaxInfoSize,offset:integer;
  253.   ch:string;
  254. begin
  255.   MaxInfoSize:=480;
  256.   SetCustom(MaxInfoSize,1,HdrSize+Width*Height);
  257.   SetImport('8-bit'); {Don't invert}
  258.   Import('');
  259.   GetRow(0,0,MaxInfoSize);
  260.   Dispose;
  261.   SetNewSize(460,100);
  262.   SetForeground(255);
  263.   SetBackground(0);
  264.   MakeNewWindow('Info');
  265.   SetCursor('Watch');
  266.   SetFont('Monaco');
  267.   SetText('With background; Left Justified');
  268.   SetFontSize(12);
  269.   MoveTo(8,10);
  270.   for i:=0 to MaxInfoSize-1 do begin
  271.     offset:=i mod 96;
  272.     if offset=0 then writeln;
  273.     ch:=chr(LineBuffer[i]);
  274.     if (offset=2) and (ord(ch)=0) then exit;
  275.     if (offset>=16) and (offset<=95) and (ord(ch)>=32) and (ord(ch)<=126)
  276.       then write(ch);
  277.   end;
  278. end;
  279.  
  280.  
  281. macro 'Import Biorad MRC 600 Z Series…';
  282. {
  283. Imports a Z series(multiple images per file) from a Biorad MRC 600
  284. confocal microscope.  The width, height and number of images are
  285. extracted from the first 3 16-bit word in the 76 byte header and
  286. the file name is extracted from bytes 18-23 of the header. This macro
  287. does not read merged pseudocolored BioRad files. Note that the Undo
  288. and Clipboard buffers must be set to 384K to work with the typical
  289. 768x512 Biorad images.
  290. }
  291. var
  292.   width,height,nImages,offset,hdrsize,i,start,picsize:integer;
  293. begin
  294.   RequiresVersion(1.50);
  295.   width:=512; 
  296.   height:=1;
  297.   offset:=0;
  298.   SetImport('8-bit'); 
  299.   SetCustom(width,height,offset);
  300.   Import(''); {Read header}
  301.   GetPicSize(width,height);
  302.   if (width<>512) or (height<>1) then begin
  303.     Dispose;
  304.     PutMessage('Please to not change width, height, etc. in the Import dialog box.');
  305.     exit;
  306.   end;
  307.   width:=GetPixel(0,0)+GetPixel(1,0)*256;
  308.   height:=GetPixel(2,0)+GetPixel(3,0)*256;
  309.   nImages:=GetPixel(4,0)+GetPixel(5,0)*256;
  310.   Dispose;
  311.   hdrsize:= 76;
  312.   picsize:=width*height;
  313.   if (width<128) or (width>2048) or (height<128) or (height>2048) or (nImages<1) or (nImages>256) then begin
  314.     PutMessage('This does not seem to be a Biorad MRC 600 Z Series file.');
  315.     exit;
  316.   end;
  317.   start:=GetNumber('Starting image:',1);
  318.   offset:=HdrSize+(start-1)*PicSize;
  319.   SetImport('8-bit; Invert'); 
  320.   SetCustom(width,height,offset,nimages);
  321.   Import('');
  322.   for i:=1 to nSlices do begin
  323.     SelectSlice(i);
  324.     Invert;
  325.     ChangeValues(0,0,1);
  326.     ChangeValues(255,255,254);
  327.   end;
  328.   ShowBioRadInfo;
  329. end;
  330.  
  331.  
  332. macro 'Import from IBAS';
  333. var
  334.    width,height,offset:integer;
  335. begin
  336.    width:=128; 
  337.    height:=1;
  338.    offset:=0;
  339.    SetImport('8-bit'); 
  340.    SetCustom(width,height,offset);
  341.    Import(''); {Read in header as an image, prompting for file name.}
  342.    width := (GetPixel(7,0)*256) + GetPixel(6,0);
  343.    height := (GetPixel(9,0)*256) + GetPixel(8,0);
  344.    Dispose(nPics);  {The ID of the last window opened = nPics.}
  345.    offset:=128;  {The IBAS offset}
  346.    SetImport('8-bit; Calibrate; Autoscale');
  347.    SetCustom(width,height,offset);
  348.    Import('');  {No prompt this time; Import remembers the name.}
  349.    Invert
  350.    SetScaling ('Bilinear');
  351.    SetScaling ('New Window');
  352.    ScaleAndRotate (0.80, 1.0, 0);
  353. end;
  354.  
  355.  
  356. macro 'Import 64x64x64x16-bit SPECT Image…';
  357. {Imports a 64x64x64x16-bit headerless SPECT image into a stack.}
  358. var
  359.   width,height,nImages,hdrsize:integer;
  360. begin
  361.   RequiresVersion(1.50);
  362.   width:=64;
  363.   height:=64;
  364.   nImages:=64;
  365.   HdrSize:= 0;
  366.   SetImport('16-bit Unsigned, Swap Bytes');
  367.   {SetImportMinMax(0,2500);} {Uncomment to fix scale}
  368.   SetCustom(width,height,HdrSize,nImages);
  369.   Import('');
  370.  end;
  371.  
  372.  
  373. macro 'Import 8-bit 3D Image…';
  374. var
  375.   width,height,offset,nImages:integer;
  376. begin
  377.   RequiresVersion(1.50);
  378.   width:=GetNumber('Width:',256);
  379.   height:=GetNumber('Height:',256);
  380.   nImages:=GetNumber('Depth(number of slices):',128);
  381.   offset:=GetNumber('Offset(header size):',0);
  382.   SetImport('8-bit'); 
  383.   SetCustom(width,height,offset,nImages);
  384.   Import('');
  385. end;
  386.  
  387.  
  388. macro 'Open with scale set to 100 pixels/mm [S]';
  389. {Example of a way to open images and have the}
  390. {spatial scale always set the same way.}  
  391. begin
  392.   Open('');           {or Import('')}
  393.   SetScale(100,'mm'); {Change as needed}
  394. end;
  395.  
  396.  
  397. macro 'Load Synergy Image';
  398. begin
  399.    SetImport('8-bit'); 
  400.    SetCustom(512,480,16384);
  401.    Import('');
  402.    ChangeValues(0,0,1);
  403.    ChangeValues(255,255,254);
  404.    SetPalette('Rainbow');
  405. end;
  406.  
  407.  
  408. macro 'Import Siemens 3D MRI…';
  409. begin
  410.   RequiresVersion(1.50);
  411.   SetImport('16-bit Signed, Swap Bytes');
  412.   {SetImportMinMax(0,3000);} {Remove comments to fix scale}
  413.   SetCustom(256,256,0,127);
  414.   Import('');
  415. end;
  416.  
  417.  
  418. procedure ImportPhotoshop(offset:integer);
  419. {
  420. Imports a 24-bit, uncompressed Photoshop
  421. TIFF file into a 3-slice stack.
  422. }
  423. var
  424.   width,height,temp,stack,i,j:integer;
  425. begin
  426.   SetImport('8-bit; Invert');
  427.   width:=GetNumber('Width:',512);
  428.   height:=GetNumber('Height:',512);
  429.   SetCustom(3*width,height,offset,1);
  430.   Import('');
  431.   temp:=PicNumber;
  432.   SetNewSize(width,height);
  433.   MakeNewStack('RGB');
  434.   AddSlice;
  435.   AddSlice;
  436.   stack:=PicNumber;
  437.   for i:=0 to width-1 do begin
  438.     for j:=1 to 3 do begin
  439.       ChoosePic(temp);
  440.       GetColumn(3*i+(j-1),0,height);
  441.       ChoosePic(stack);
  442.       ChooseSlice(j);
  443.       PutColumn(i,0,height);
  444.     end;
  445.   end;
  446.   SelectSlice(1);
  447.   SelectPic(temp);
  448.   Dispose;
  449. end;
  450.  
  451. macro 'Import Photoshop 2.0 24-bit TIFF';
  452. begin
  453.   ImportPhotoshop(196);
  454. end;
  455.  
  456. macro 'Import Photoshop 2.5 24-bit TIFF';
  457. begin
  458.   ImportPhotoshop(598);
  459. end;
  460.  
  461. macro 'Import LUT…';
  462. {Imports a 256 x 3 x 8-bit look-up table located 'offset'
  463.  bytes from the beginning of a file. Use an offset of 32 for
  464.  LUTs created using Image's Save As command and an
  465.  offset of 0 for Exported LUTs.}
  466. }  
  467. var
  468.   offset,i:integer;
  469. begin
  470.   offset:=0; {Use 32 for LUTs created using Save As}
  471.   SetImport('8-bit'); 
  472.   SetCustom(256,3,offset);
  473.   Import(''); {Read LUT as an image}
  474.   for i:=0 to 255 do begin
  475.     RedLut[i]:=GetPixel(i,0);
  476.     GreenLut[i]:=GetPixel(i,1);
  477.     BlueLut[i]:=GetPixel(i,2);
  478.   end;
  479.   UpdateLUT;
  480. end;
  481.  
  482. macro 'Open List of Files…';
  483. var
  484.   width,height:integer;
  485.   name,ch:string;
  486. begin
  487.   width:=256;
  488.   height:=3;
  489.   SetCustom(width,height,0);
  490.   Import('File List'); {Read text as an image}
  491.   repeat
  492.     name:='';
  493.     repeat
  494.        ch:=chr(GetPixel(CharLoc mod width,CharLoc div width));
  495.        if ch>=chr(32) then begin
  496.          name:=concat(name,ch);
  497.        end;
  498.        CharLoc:=CharLoc+1;
  499.      until (ch=13) or (CharLoc>=width*height);
  500.      ShowMessage(name);
  501.      wait(1);
  502.   until (name='') or (CharLoc>=width*height);
  503. end;
  504.  
  505. macro 'File processing test';
  506. begin
  507.   Open('hd400:images:A');
  508.   Invert;
  509.   Save;
  510.   Close;
  511.   Open('hd400:images:B');
  512.   Invert;
  513.   Save;
  514.   Close;
  515.   Open('hd400:images:C');
  516.   Invert;
  517.   Save;
  518.   Close;
  519. end;
  520.  
  521.  
  522.  
  523.  
  524.